Inelastic X-ray scattering study of the collective dynamics in liquid sodium. 
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Inelastic X-ray scattering data have been collected for liquid sodium at T = 390 K, i.e. slightly above 
the melting point. Owing to the very high instrumental resolution, pushed up to 1.5 meV, it has been 
possible to determine accurately the dynamic structure factor, S(Q,lj), in a wide wavevector range, 
1.5 15 nm~^, and to investigate on the dynamical processes underlying the collective dynamics. A 
detailed analysis of the lineshape of S(Q,uj), similarly to other liquid metals, reveals the co-existence 
of two different relaxation processes with slow and fast characteristic timescales respectively. The 
present data lead to the conclusion that: i) the picture of the relaxation mechanism based on a 
simple viscoelastic model fails; ii) although the comparison with other liquid metals reveals similar 
behavior, the data do not exhibit an exact scaling law as the principle of corresponding state would 
predict. 



I. INTRODUCTION 

Liquid metals are well known to exhibit remarkably 
pronounced inelastic features (Brillouin peaks) in their 
density fluctuations spectra up to very high Q values 
(well beyond the hydrodynamic regime). For this main 
reason, in the last decades, the dynamic structure factor 
of these systems has been used as a benchmark for the 
comprehension of the mechanisms underlying the atomic 
motions at the microscopic level in liquids. In particu- 
lar, in the special case of alkali metals, collective oscilla- 
tory modes persist down to wavelengths of the order of 
one or two interparticle distances, a feature which makes 
them ideal candidates to test different models for collec- 
tive properties at finite wavevectors. 

In the last thirty years many experimental investiga- 
tions have been performed on these systems by means 
of inelastic neutron scattering (INS), which was, up to a 
few years ago, the only tool adequate to investigate the 
dynamic of condensed matter in the mesoscopic wavevec- 
tor region. Indeed, since the very first inelastic neutron 
scattering experiments by Copley and Rowe in liquid 
rubidium, many similar studies have been reported for 
similar systems: lithium sodium j^, potassium [Q, 
cesium 1^ and again rubidium . In the most favorable 
cases (for example in cesium), important features such 
as the so called positive dispersion of sound, namely an 
upwards deviation from the low-Q, linear dispersion be- 
havior, have been reported. Moreover, a description of 
the coherent dynamic structure factor has been given in 
terms of a single relaxation process (viscoelastic model) . 
On the other hand, in some of the above mentioned sys- 
tems, an accurate investigation of the collective prop- 
erties is prevented by two drawbacks of INS technique, 
namely the presence of both coherent and incoherent con- 
tribution to the inelastic scattering cross section, and 
the occurrence of kinematic restrictions which confine the 



minimum accessible wavevector above a threshold value 
which is already outside the linear dispersion regime. 
This is particularly the case of lighter elements such as 
lithium and sodium. 

Despite the lack of experimental information on the 
collective properties in these systems, significant ad- 
vances have been done by means of computer simulation 
techniques, where the afore mentioned restrictions do not 
apply. Consequently, many numerical studies have been 
reported aiming at the study of the density fluctuation 
spectra of liquid metals ||7|-p^. Moreover, by these nu- 
merical techniques one can access a wider set of correla- 
tion functions, while an inelastic scattering experiment 
basically probes the density autocorrelation function. 

On the theoretical side, thanks to development of ap- 
proaches such as the memory function formalism, the re- 
laxation concept, the kinetic theory, it has been 
built up a framework able to account for the behavior 
of the afore mentioned correlation functions. Along this 
line, a general picture of the dynamics of the density fluc- 
tuations in simple liquids gradually emerged. In particu- 
lar it is now established that the decay of the density au- 
tocorrelation function occurs through mechanisms char- 
acterized by different time scales. In addition to a first 
mechanism related to the coupling of the density and the 
temperature modes (thermal relaxation), one must allow 
for the dynamical effects of the stress correlation func- 
tion (viscous relaxation). This latter process is thought 
to proceed throughout two distinct relaxation channels, 
active over quite different time scales: a first rapid de- 
cay, customarily ascribed to the interaction of each atom 
and the "cage" of its neighboring atoms, is followed by 
a slower process which yields a long lasting tail. This 
mechanism stems from slow, temperature-dependent re- 
arrangements which eventually, in those systems capable 
to sustain supercooling, may cause a structural arrest 
(glass transition) . This structural relaxation process has 
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widely been studied within the mode coupling formal- 
ism. Recently, in the specific case of liquid alkali metals, 
remarkable attempts have been made to set up self con- 
sistent approaches which include all these mechanisms 
from the very start ||2|,|2|l. 

This theoretical framework has been tested on a num- 
ber of numerical studies (see, for example, [Q). On the 
experimental side, however, the lineshapes extracted by 
INS in those systems where one may access a significant 
{Q — E) region did not allow to discriminate between dif- 
ferent models including all the relevant relaxation pro- 
cesses 1^. 

Although trough of INS enormous advances in the com- 
prehension of the collective properties of condensed mat- 
ter have been reported, some basic aspects remain still 
unsettled. In particular, due to the previously mentioned 
features of INS, sufficiently accurate measurements of the 
coherent dynamic structure factor, S{Q,uj) have hardly 
been reported, with the result that a quantitative de- 
scription of the relaxation dynamics in monoatomic fluids 
was mainly based on simulation data. 

Only in the recent past, the development of new syn- 
chrotron radiation facilities opened the possibility of us- 
ing X-rays to measure S{Q,lu) in the non-hydrodynamic 
region; in this case the photon speed is obviously much 
larger than the velocity of the excitations, and no kine- 
matic restriction occurs. Moreover, in a monatomic sys- 
tem, the X-ray scattering cross section is purely coher- 
ent, and consequently can be directly associated with the 
coherent dynamic structure factor |l4|j2^ ]. After some 
early, pioneering data ||2^], IXS rapidly developed in the 
last decade, becoming fully operational with the advent 
of the third generation sources 1 28 2^ . The applicability 
of IXS to the study of collective and single particle dy- 
namics of light liquid metals has been then exploited in 
several works performed at the beamline ID16 of the eu- 
ropean synchrotron radiation facility (ESRF) in Greno- 
ble [^0|-p2| , thus showing the complementarity of IXS and 
INS. More recently, IXS has been pushed to the top of 
its theoretical performance in terms of flux and resolu- 
tion, and therefore new and more accurate IXS studies on 
liquid metals allowed to explore the relaxation processes 
underlying the microscopic dynamics in monoatomic flu- 
ids [^3|-^. In particular, it has been possible to detect 
experimentally and to quantify the presence of the re- 
laxation scenario predicted by the theory of simple liq- 
uids. In detail, the data have been analyzed within a 
generalized hydrodynamic approach ||lj,|l^ and, follow- 
ing a memory function formalism | [T^ in the mesoscopic 
wavevector region, it has been possible to detect the pres- 
ence of two viscous relaxation processes in addition to the 
thermal process [|4j. These results, obtained with a phe- 
nomcnological ansatz for the memory function originally 
proposed in a theoretical work by Levesque et al. [ pT[ , 
has then been deeply discussed in order to clarify the 
physical origin of all these dynamical processes . 



In this work, we report the results of a very high res- 
olution determination of the dynamic structure factor in 
liquid sodium at the melting point at several fixed Q val- 
ues. Several approximation are adopted for the analysis 
of the experimental data, and the number and the role 
of the involved relaxation processes is discussed in detail. 
Finally, a comparison with similar experimental data on 
different liquid metals is presented. 



II. THE EXPERIMENT 

In liquid sodium, the coherent and incoherent INS 
cross sections are almost equivalent, and the adiabatic 
speed of sound is about 2500 m/s. As a result, a study 
of the collective properties in the Q-region characterized 
by a linear dispersion (the first Brillouin pseudo-zone) 
is quite hard to perform ||^. The use of IXS to these 
systems is, therefore, particularly advantageous and, as 
shown in ref. , the inelastic features of the dynamic 
structure factor can easily be detected even around the 
Q position of the first sharp diffraction peak {Qm ~ 20 
nm""'^). 

Data have been collected for liquid sodium at the new 
IXS beamline ID28 of ESRF (which is now fully oper- 
ative), allowing an energy resolution of AE = 1.5 meV 
FWHM. The experiment has been performed at fixed 
exchanged wavevector over a Q— region 1.5 15 nm~^. 
Each energy scan (—50 < E < 50 meV) took about 300 
minutes, and has been repeated for a total integration 
time of about 500 seconds/point. A five analyzers bench, 
operating in horizontal scattering geometry, allowed us to 
collect simultaneously photons at five different values of 
exchanged wavevector Q for each single scan. As sample 
environment we used a setup similar to the one adopted 
in previous experiments performed on liquid lithium , 
namely an austenitic steel cell heated by thermal con- 
tact with a resistor connected to a voltage regulated sup- 
ply. The sample length has been optimized by matching 
the absorption length of sodium at the working energy 
21747 eV (corresponding to the Si (11 11 11) reflection 
in backscattering geometry) which is about 5 mm. 

In Fig. 0, a selection of experimental spectra is re- 
ported. Since the incident fiux on the sample varies with 
time, data have been normalized to the monitor. More- 
over, since the five analyzers have different scattering effi- 
ciencies, to extract S{Q, uj) from the raw spectra the data 
have been put in absolute units using the knowledge of 
the lower order frequency momenta of the dynamic struc- 
ture factor: 
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where the tilde indicates the true quantum dynamic 
structure factor, to distinguish it from its classical repre- 
sentation. 
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FIG. 1. Selected IXS energy spectra of liquid sodium at 
T = 390 K and fixed values of exchanged wavevector (open 
dots o). The data have been normalized utilizing the sum 
rules according to the procedure described in the text. 

This procedure can be summarized as follows. The 
actual experimental spectrum I{Q,ijj) is proportional to 
the convolution of the scattering law S{Q,uj) with the 
resolution R{lo): 



I{Q,uj)^E{Q) / duj'S{Q,uj')R{uj -uj') 



(3) 



In terms of the first two spectral moments of / and R: 
^fl^ ^B^^ easily deduces that 
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Consequently, the normalized spectrum (still affected by 
the resolution broadening) reads: 
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The reliability of this procedure will be discussed in the 
next section. 



III. DATA ANALYSIS 

The data analysis has been performed following a 
memory function approach, similar to the one adopted 
for lithium and aluminum data p3|-^. 

Within the generalized Langevin equation formalism, 
it is possible to express the classical dynamic structure 
factor in terms of a complex function, M(Q,t), (memory 
function) related to the interaction details. In particular, 
S{Q,ll!) can be expressed in terms of the real and imag- 
inary part of M{Q,uj), the Fourier-Laplace transform of 
M{Q,t) as [0: 



S{Q,u;) 
S{Q) 



Tr-^Lui{Q)M'{Q,uj) 



uj^ ~ujI{Q)+luM"{Q,uj) + LuM'iQ.Lo) 



(6) 

where the quantity t^oiQ) — mS{Q) related to the 
generalized isothermal sound speed through the relation 

ctiQ)=MQ)/Q- 

In order to use the above expression to reproduce the 
experimental spectra, I^iQ,!^), it has to be modified to 
be quantum-compliant with the detailed balance condi- 
tion, and finally must be convoluted with the instrumen- 
tal resolution R(uj). Utilizing one of the most common 
quantum transformation, one finds that 



huj'/KT 

2 _ ^-huj' /KT 



SiQ,Lu')R{LU-Lu')dLj' (7) 



The above expression can finally be used as a fitting 
function to the experimental data to obtain the rele- 
vant relaxation parameters, i.e. relaxation times and 
strengths. In this procedure, the values of S{Q) arc ex- 
tracted from the experimental data exploiting Eq. (^. 
As memory function, M{Q,t), we utilized a multiple ex- 
ponential ansatz. In Ref. js^, indeed, it has been shown 
that in a similar system, namely liquid lithium: i) at 
least three distinct relaxation processes (one thermal and 
two viscous) are necessary to reproduce the IXS line- 
shape; and ii) a multi-exponential shape is an acceptable 
approximation of the actual memory function. Conse- 
quently, the total memory function reads 

Af(Q,i)-(7-l)c.2(Q)e-^-'5'* (8) 
+ A2(Q) U(Q)e-*/"°(^) + (1 - A(Q))e-*/""(^)' 

The first term of the above equation comes from the cou- 
pling between thermal and density degrees of freedom, 
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and depends on the thermal conductivity Dt and the 
specific heat ratio 7. The remaining two contributions in 
Eq. are the genuine viscous processes. Through the 
fitting procedure, the unknown parameters r^, Tq, (the re- 
laxation times), A{Q) (the relative weight of the viscous 
processes), and A^(Q) (the total viscous strength), can 
be deduced. 
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FIG. 2. IXS spectra (open dots o) fitted with one (thin 
line) or two (thick line) relaxation process. Experimental res- 
olution is also indicated (dotted line) 

In Fig. 1^, we report the outcome of this fitting proce- 
dure, i.e. a comparison between Ij^{Q,uj) and In{Q,(^)- 
Both the results of a one and a two time ansatz are shown. 
As in the case of lithium and aluminum, a single re- 
laxation time description (viscoelastic model) appears to 
yield a rather poor reproduction of the data. In particu- 
lar, at low Q's the simultaneous presence of two processes 
is crucial even to catch just the qualitative features of the 
spectra. 

As previously mentioned, once that the experimental 
values of S{Q) are introduced in Eq. the only free fit- 
ting parameters are the relaxation times and strength. A 



consistency check has however been obtained comparing 
these values of S{Q) (sum rules) with the one calculated 
from ujq{Q) left as free parameter in Eq. (^. 



0.30 



0.25 



0.20 



5? 



0.15 



0.10 



0.05 



0.00 




6 8 10 12 14 16 

Q (nm ') 

FIG. 3. Values of S(Q) from the spectral moments of the 
IXS signal (full line), compared with the values taken from 
the fit (full circles o) and from literature (dotted line • • ■) 



0.06 



0.04 - 



0.02 
1.0 

0.8 

0.6 

0.4 

0.2 

0.0 



i • • • 



10 12 14 



16 



Q(nm ) 



FIG. 4. Main frames: Fast and slow relaxation times as 
obtained by the fitting procedure. Insets: product of the 
relaxation times with the position of the inelastic peak of the 
classical density correlation function, as deduced by the fit. 

In Fig. H the results are shown, together with the ex- 
perimental determination of S{Q) from neutron diffrac- 
tion data pq l . 

In Fig. W we show the two relaxation times, Tq and 
as obtained by the fitting procedure. The slower process 
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is almost constant, with an abrupt increase below Q « 6 
nm~^. Such an increase is probably an artifact, since 
it occurs as the timescale approaches the inverse of the 
resolution width, here reported as a dashed line. Conse- 
quently the determination of the a relaxation time below 
this Q value is no longer reliable. The fj, relaxation pro- 
cess shows instead a slightly decreasing Q dependence. 
In the insets, we report lut as function of Q for the two 
processes, where omega is the dominant frequency of the 
density fluctuations at a given Q (i.e. the frequency of 
the Brillouin component). The quantity lot gives im- 
portant indications about the role of the two processes. 
The faster one, indeed, has values of lot always smaller 
than one, so that it is mainly responsible for the acous- 
tic damping (Brillouin linewidth). Moreover, u}[Q)t^{Q) 
shows an increase in the same wavevector region where 
the sound dispersion occurs (approximately, up to Q val- 
ues of one half of the position of the first sharp diffraction 
peak), so that an increase of the apparent peak position 
ruled by the strength of the fast process has to be ex- 
pected. The Lo{Q)Ta{Q) is always larger than one. Con- 
sequently the alpha process mainly control the sharper 
portion of the quasielastic scattering, while it only gives 
a moderate contribution to the speed of sound over all 
the explored Q range (the strength of the a process is 
negligible with respect to that of the fast relaxation). 
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FIG. 5. Inset: dispersion curve uJi{Q) from the deconvo- 
luted lineshapes (full dots •) compared with the low fre- 
quency, isothermal, behavior (dotted line ■ ■ •)• Main plot: 
speed of sound in the same notation of the inset. The in- 
stantaneous value Coo(Q) (full line) is also reported. The 
co{Q — > 0) limit of the isothermal value, from compressibility 
data, is shown as dotted arrow. 

The above mentioned increase of the speed of sound 
is illustrated in Fig. ^, where coiiQ), the maximum of 
the longitudinal current correlation spectrum J^{Q^ lu) = 
uj'^ /Q^S{Q, uj) is reported in the main frame, and the as- 
sociated sound speed uJi{Q)/Q is reported in the inset. 



The positive dispersion turns out to be of the order of 
20%, similarly to what is found in other liquid metals 
Q. In contrast with the results reported in Ref. [ p2[ , 
we find that the isothermal value of the sound speed is 
not yet reached at the minimum wavevector probed in 
our experiment, i.e. 1.5 nm~^. At this low Q resolution 
effects are particularly important: even with the present 
resolution (1.5 meV FWHM) the Brillouin components 
of the S{Q,oj) appear as shoulders of the quasielastic 
signal, therefore we believe that this discrepancy may 
be ascribed to the lower resolutions of the previous ex- 
periment that prevents an accurate determination of the 
actual lineshape at such small Q value. Both the isother- 
mal Ct{Q) (low frequency) ]3^ ] and the infinite frequency 
CooiQ) values of the sound speed have been reported as 
obtained by the fit. This latter quantity is always higher 
than the maximum value reached by the apparent sound 
speed, consistently with the behavior of the relaxation 
times: as shown in FigjH the faster process (which is the 
one with the dominant strength) has a relaxation time 
which never exceeds the dominant timescale of the den- 
sity fluctuation at the probed wavevector. This notwith- 
standing, the actual values of Coo(Q) reported in Fig. |^ 
can partially be affected by the oversimplifications given 
by the assumption of exponential lineshape for the mem- 
ory functions fs^ . 

We conclude this investigation with the comparison of 
the S{Q,Lu) of three different liquid metals. In Fig. ^ 
the IXS data on liquid sodium arc compared with those 
of lithium and aluminium by introducing a set of length, 
mass and time units [I* ,m* ,t*] for each system. Such 
a procedure can be considered meaningful as far as is 
valid a corresponding state principle. As length unit, I*, 
we have chosen /* = Q^^, the inverse of the position of 
the static structure factor S(Q). The time unit, t* , has 
been instead chosen as t* = ^ I -^l* , with Tm the melt- 



ing temperature and m the atomic mass. In the upper 
panel the raw IXS signal, scaled on this units is reported. 
The Q values have been selected among the available 
ones as closer as possible for three different regions of 
the dispersion curve. It is worth to point out that all the 
spectra are resolution broadened and, moreover, that the 
three experiments have been performed in different con- 
dition of resolutions. To bypass the instrumental resolu- 
tion bias, one can compare the lineshape obtained by the 
fits, i.e. the classical deconvoluted dynamical structure 
factor that better reproduces the experimental data after 
the transformation of Eq. (^. The result of such a com- 
parison is reported in the lower panel of Fig. |^. The alu- 
minum spectra clearly have a different lineshape, while 
lithium and sodium are much more similar when scaled 
in their respective absolute units. This notwithstanding, 
the lineshapes of the two alkali systems do not coincide, 
implying that the principle of corresponding state is not 
entirely valid for these systems. One comes to the same 
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conclusion even looking at some basic quantities such as 
the sound speed: the scaling factors between lithium and 
sodium, for example, are not the same if one looks at 
Ct{Q 0) or Coc(Q 0). Also the static structure fac- 
tor S{Q), rather of being universal quantity, depends on 
the considered system. 



terizing both alkali and non- alkali simple liquids, some 
important quantitative differences are found: in particu- 
lar, a scaling law among alkali systems seems to be unap- 
plicablc, in contrast to some previous findings obtained 
by computer simulation data. 
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FIG. 6. Comparison of the IXS spectra of sodium, lithium 
and aluminium reported in scaled units, as defined in the 
text. Best fitted lineshape (two time relaxation model) are 
also shown, b) same as a) but with a classical, resolution 
deconvoluted, lineshape as obtained from the values of the 
parameters deduced by the fitting procedure. 



IV. CONCLUSIONS 

In conclusion, we presented a very high resolution ex- 
perimental study of the collective dynamic structure fac- 
tor in liquid sodium at the melting temperature. In order 
to account for all the spectral details, beside the usual 
thermal process coupling temperature and density, at 
least two additional viscous relaxation processes have to 
be invoked. The viscoelastic approximation do not give, 
indeed, a satisfactory description of our data. These re- 
sults have finally been compared with those obtained in 
other liquid metals. Although the mechanisms ruling the 
dynamics turn out to lie in a general framework charac- 
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